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ABSTRACT 

We present gravitational lens models for 20 strong gravitational lens systems observed 
as part of the Sloan WFC Edge-on Late-type Lens Survey (SWELLS) project. Fifteen 
of the lenses are taken from paper I while five are newly discovered systems. The 
systems are galaxy-galaxy lenses where the foreground deflector has an inclined disc, 
with a wide range of morphological types, from late- type spiral to lenticular. For each 
system, we compare the total mass inside the critical curve inferred from gravitational 
lens modelling to the stellar mass inferred from stellar population synthesis (SPS) 
models, computing the stellar mass fraction /* = Msps/Mi ens . We find that, for the 
lower mass SWELLS systems, adoption of a Salpeter stellar initial mass function (IMF) 
leads to estimates of /* that exceed 1. This is unphysical, and provides strong evidence 
against the Salpeter IMF being valid for these systems. Taking the lower mass end of 
the SWELLS sample (asm < 230 km s -1 ), we find that the IMF is lighter (in terms 
of stellar mass-to-light ratio) than Salpeter with 98% probability, and consistent with 
the Chabrier IMF and IMFs between the two. This result is consistent with previous 
studies of spiral galaxies based on independent techniques. In combination with recent 
studies of massive early-type galaxies that have favoured a heavier Salpeter-like IMF, 
this result strengthens the evidence against a universal stellar IMF. 

Key words: galaxies: spiral - galaxies: fundamental parameters - stars: mass function 
- gravitational lensing 



1 INTRODUCTION 

The stellar initial mass function (IMF) is a fundamental 
property of a simple stellar population, and plays a key role 
in many astrophysical problems. For instance, the determi- 
nation of galaxy stellar masses and star formation rates are 
uncertain at the factor of ^ 2 level due to uncertainty in th e 
IMF fe.g. lConrov. Gunn. fe Wbiitel 20091 ; Fanter et alJ2007h . 
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Observations of the solar neighbou rhood favor an IMF 
with a power law slope at h igh mass es jSaipeterlll955h. and 
a turnover at low masses (iKroupal l200ll : Chabrier 120031 ). 
While other spir al galaxies have IMFs consistent with that of 
the Milky Way (iBell fe de JondlioOll: iBershadv et alJl201ll : 



button et all l2011al lbl: ISuvu et al.1 1201 lh . there is increas- 
ing evidence from a variety of observations that the IMF 
is non-universal. Specifically, massive early-type galaxies at 
low redshifts r equire IMFs that are heavier than in the 
Milky Way fe.g.lSpiniello et al.ll201ll : ISonnenfeld et alJ201ll : 
lAueer et alj|2010bh . and po ssibly heavier than Salpeter i n 
the most massive ellipticals flvan Dokkum fe Conrovll2010h . 
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Other lensing and dynamical studies of e arly-type galaxies 
also suggest that th e IMF may be heavy ([Treu et alj|20ld : 
iThomas et al.|[201lh although these observations may also 
be explained by a trend in the properties of the dark matter 
haloes. Recently, variations in the IMF of early- t ype g alax- 
ies have also been reported bv | Cappellari et all (|2012h and 
button. Mendel, fe Simardl (|2012h . 

One of the simplest ways of constraining the galaxy 
averaged IMF is through comparing total masses (within 
some aperture), as derived by kinematics and/or strong 
gravitational lensing and referred to as gravitational masses 
henceforth, with stellar masses de rived from stella r pop- 
ulation synthesis (SPS) models (iCappellari et all |2006; 
iFerreras et al.| [201Q; Th omas et al.l l201lh . While this method 
typically only provides strong upper limits to the stellar 
mass due to the unknown dark matter fraction, as well as 
the unknown mass in gas, dust etc., the upper limits are 
robust. 

Dark matter fractions in galaxies are expected to in- 
crease with radius, so the strongest constraints on the IMF 
(from total mass measurements) come from measuring mass- 
to-light ratios at small galactic radii. However, at small radii 
kinematics are harder to interpret, as there is usually a mix 
of rotation and dispersion. This is where strong gravitational 
lensing is particularly useful and unique, because it mea- 
sures projected mass independent of the dynamical state of 
the deflector. 

Assuming a universal IMF, the dark matter frac- 
tions within the effective radius of early-type galax- 
ies decrease with decreasing mass (down to a stel- 



lar mass of ^ 10 Mq) (e.g. [Padmanabha n et al 



Gallazzi et al .1120061: lAuger et al.ll2010al : button et al.1 



Barnabe et al 



2004 



2011g 



1201 ll ). Thus lower mass elliptical galaxies are 



expected to place stronger constraints on the IMF than 
higher mass elliptical galaxies (although high mass systems 
can place strong con straints as well through dynamics, e.g. 
ISpiniello et al.l 1201 lh . If the bulges of spiral galaxies fol- 
low the same scaling relations as elliptical galaxies then 
we would expect that bulges of mass ~ 1O 1O M - where 
the dark matter fraction is minimised - would provide the 
strongest constraints on the IMF if the IMF is univers al (e.g. 
iBastian. Covev. fc Meverl l2010l : iKrouoa et all l2Qllh . Con- 
versely, if the IMF is not universal, measuring it as a function 
of galaxy mass, morphology, and environment should help 
understand the physical mechanisms that determine it. 

Measuring the absolute normalisation of stellar IMF as 
a function of galaxy type and mass is one of the main goals of 
the Sloan WFC Edge-on Late- type Lens Survey (SWELLS; 
Treu et al. 2011, Dutton et al. 2011b; hereafter papers I and 
II, respectively). 

In this paper, we present Hubble Space Telescope (HST) 
multicolour imaging of the complete SWELLS sample, in- 
cluding new data from cycle- 18 observations, and use strong 
lensing models to constrain the stellar initial mass function 
of spiral galaxies. In particular we use precision measure- 
ments of the total mass within the critical curve and there- 
fore place an upper limit on the amount of stellar mass. 
We combine this inferences with stellar masses derived from 
stellar population synthesis models to obtain strong upper 
limits on the normalisation of the stellar IMF. 

The galaxies in the SWELLS sample have total stellar 
masses (assuming a Chabrier IMF) of M* ~ 4-3OxlO lo M , 



and bulge masses of ~ 2 - 10 x 1O 1O M (jTreu et al.ll201lh . 
In order to explore the variation of the IMF with stellar 
mass, we combine our sample with the more more massive 
ellipt ical galaxy lenses from the SLACS survey (| Auger et al.l 
l20Q9h with M* ~ 10 - 50 x 1O 1O M . 

Throughout this paper, we assume a flat ACDM cosmol- 
ogy with present day matter density, Q m = 0.3, and Hubble 
parameter, Hq — 70kms _1 Mpc -1 . 



2 THE DATA 

In this section we briefly summarise the sample selection 
procedure (§ I2.H more details can be found in paper I) 
and then present the new cycle- 18 H ST observations in Sec- 
tion E21 



2.1 Sample Selection 

The sample of spiral lens galaxy candidates was selected 
from the SDSS database as discussed in paper I. Briefly, 
the spectroscopic database was searched for composite spec- 
tra, consisting of foreground spectrum superimposed with 
multiple emission lines at a higher red shifts, following the 
proce dures developed for the SLACS (iBolton et al.l 20061. 
20081 . hereafter B06 and B08) and BELLS (jBrownstein etaD 
20121 ) surveys. In addition, however, in order to identify in- 
clined late-type galaxies suitable for rotation curve measure- 
ments, the deflector galaxies were selected to have axis ratios 
q = b/a < 0.6 with a preference for galaxies with edge-on 
discs and emission lines indicating star formation. A priori 
lensing probabilities were estimated based on the SDSS stel- 
lar velocity dispersion and stellar Einstein Radius (paper I) . 



2.2 HST observations and surface photometry 

The highest probability lens candidates were followed up 
with HST in supplemental cycle 16s program (GO- 11978; PI: 
Treu) and are presented in paper I. In addition, during the 
same program, multiband imaging was obtained for disc lens 
galaxies identified as part of the SLACS survey (GOs 10174, 
10494, 10798, 11202 PI: Koopmans; GO 10587, 10886, PI: 
Bolton). Lower probability systems were included in a cy- 
cle 18 program (GO-12292; PI: Treu), aiming to extend the 
sample to later type spirals, and are presented here for the 
first time. Two candidates confirmed from the ground with 
adaptive optics (paper I) were also included in the cycle 18 
program to collect optical multiband photometry. Twenty 
one candidates were observed in cycle- 18. Eighteen of those 
were observed with the refurbished Advanced Camera for 
Surveys (ACS) for one orbit each, split equally between fil- 
ters F435W and F814W. Two dithered exposures were ob- 
tained for each filter to aid in cosmic ray and defect removal 
and to help improve the sampling of the PSF. Four candi- 
dates, including one observed with ACS were observed with 
the infrared channel of the Wide Field Camera 3 for a full 
orbit through filter F160W. Each orbit was split into four 
dithered exposures to aid in cosmic ray and defect removal 
and help improve the sampling of the PSF. The HST images 
were reduced using standard multidrizzle techniques. As ex- 
pected the fraction of grade 'A' and 'B' lenses was lower for 
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Figure 1. Images of the 21 candidate lens systems oberved with the HST during Cycle 18; the grade of each system is provided in 
parentheses next to the name. 
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cycle 18 than for previous cycles, 6/21 and 10/21 respec- 
tively; images for these 21 systems are displayed in Figure 

m 

In total, after HST imaging, the SWELLS sample con- 
sists of 20 secure lenses (grade 'A'), 6 probable lenses (grade 
'B'), and 12 possible lenses (grade 'C'). Grade 'A' lenses 
are defined as having unambiguous multiple images success- 
fully reproduced by a reasonably simple lens model. Grade 
'B' lenses are defined as having evidence for strong lensing, 
even though the quality of the data is not sufficient to iden- 
tify unambiguously multiple images and/or construct a lens 
model. This definition is quite conservative and we expect 
that with deeper data or narrow band imaging data to facili- 
tate removal of the foreground deflector light most 'B' grade 
systems would be confirmed as lenses. Grade 'C' are possible 
lenses where we expect that even with significantly better 
data only a minority of the systems would be confirmed as 
lenses. A colour montage of the 20 grade 'A' lenses from 
SWELLS is shown in Figure El 

For this paper we focus our analysis on the grade 'A' 
lenses, plus one grade 'B' lens (J1422+4134), for which we 
were able to obtain a plausible lens model. However, in light 
of the uncertain identification of this one system, we will 
use only the grade 'A' lenses for our inference on the stellar 
initial mass function. We also excluded one grade 'A' lens, 
J1037+3517, since the morphology of the lens is complex 
(involving a merger) and the modelling approach of this pa- 
per was unlikely to be successful for this system. A summary 
of the basic properties of the SWELLS lenses together with 
the inferred parameters of the gravitational lens models is 
given in Table [TJ A summary of the basic properties of the 
SWELLS HST Cycle 18 targets is given in Table H 

Multiband surface photometry was derived as described 
in paper I, by fittin g two-component (b ulge+disc) models 
using SPASMOID (jBennert et al.l l201lh . The model con- 
sists of an n — 4 Sersic profile (i. e. a de Vaucouleurs profile, 
ISersidll968l :lde~ Vaucouleursll 1 948h for the bulge and an n — 1 
exponential disc component. Estimates of the model param- 
eters, including stellar masses and mass-to- light ratios, are 
given in Tables G3 [H and \5\ 



3 GRAVITATIONAL LENS MODELS 

In this section we describe the procedure used to ob- 
tain gravitational lens models for each system. We con- 
sider relatively simple mass models for the deflector, de- 
scribe d as a singular isothermal ellipsoid plu s external 
shear ([Kormann, Schneider, fc Bartelmannll 1994 ) . This sim- 
ple model is known to provide a good descriptio n of (at least 
early-type) galaxy-scale gravitational len ses (|Tre 3 12010I : 
iKoopmans et al.ll2009l ; lBarnabe et al.ll2009h , and should pro- 
vide sufficiently robust estimates of the physical quantities 
relevant for this present analysis. That is, the mass enclosed 
within the critical curve, and the position angle and the el- 
lipticity of the m ass distribution. The flat ro tation curves 
of spiral galaxies flvan Albada fe Sancisilll986h suggest that 
this model is appropriate for spiral lenses, at least to first 
order. 

The modelling of the SWELLS lenses is complicated 
compared to typical gravitational lens modelling, for several 
reasons. The first challenge is the disentanglement of the 



light of the source and that of the detector. Our strategy to 
address those is discussed in § 13.11 An additional challenge 
is the proper estimation of model parameter uncertainties. 
This is discussed in § 13.21 where we present our mass mod- 
els and inference strategy. Finally, in § 13.31 we discuss each 
individual system listing any additional complexities. 

3.1 Disentangling the source and deflector light 

The first step is the selection of the band which contains 
the highest fidelity image of the lensed features, and will 
therefore be used to produce the lens model. In most cases 
this is the bluest image available, as the lensed features tend 
to be brightest at the blue end of the spectrum. In the cases 
where the lens galaxy contains significant amounts of dust, 
a redder band was selected to minimise complications from 
the dust. 

Most significantly, the small Einstein radii of the lenses 
cause the lensed features to lie in regions where the surface 
brightness of the lens galaxy is non-negligible. This makes 
lens galaxy subtraction non-trivial. In principle, we should 
fit the light profile of the lens galaxy simultaneously with 
any lensing parameters and source surface brightness profile 
parameters. We have attempted this, but found that it was 
too computationally intensive to be practical. 

In addition, the light profiles of the lens galaxies cannot 
be desc ribed adequa tely by simple models such as the Sersic 
profile (|Sersidll968h due to the presence of dust lanes, star 
forming regions, and other complexities. Attempting to use 
such a profile would result in lens galaxy subtractions with 
significant residuals that would interfere with the modelling 
of the lensed features. 

3.1.1 Rotation Subtraction for Masking 

To proceed with the lens galaxy subtractions, we adopted 
the following procedure. Firstly, for exploratory purposes, 
we take advantage of the near- symmetry of the lens galaxy: 
by subtracting the image from a 180-degree rotated version 
of itself, the lensed features, often asymmetric, are revealed. 
The centre point around which the rotation takes place can 
be left as a free parameter, and optimised to find the centre 
of rotation which makes the rotated image best match the 
original image. 

These rotation-subtracted images enable us to get an 
initial estimate for the position and morphology of the lensed 
features, particularly the identification of any faint counter- 
images. However, these images cannot be used as data for 
modelling because they contain both positive and negative 
reproductions of the arc. 

In practice, we used the rotation-subtracted images to 
inform the placement of a masked region, covering the lensed 
features. With these features masked out, we fit a complex 
many-parameter model for the surface brightness profile of 
the lens galaxy, and subtract that model from the data, giv- 
ing the final data for lens modelling. 

3.1.2 Multi- Gaussian Subtraction 

With the mask in place, we fit the remaining pixels using 
the code of iBrewer et al.l (|2011ah , with all Sersic indices set 
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Figure 2. Images of the 20 grade 'A' lens systems that comprise the SWELLS sample. 
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Table 1. Summary of basic properties of modelled grade 'A' lenses (plus the modelled grade 'B' lens J1422+4134) from the SWELLS 
sample of spiral lens galaxies. Col. 1 lists the lens ID; Cols. 2 and 3 the coordinates (J2000 in degrees); Cols 4 and 5 give deflector and 
source redshifts; Col. 6 the Einstein Radius of the lens model; Col. 7 the velocity dispersion of the Singular Isothermal Ellipsoid plus 
external shear (SIE+7) lens model; Col 8. the axis ratio of the lens model; Col. 9 the projected mass within the critical curve of the SIE 
lens model; Cols. 10 Sz 11 the ratio between the stellar mass and lensing mass within t he critical curve, a ssuming Chabrier and Salpeter 
IMF, respectively; Col. 12 is the reference for the lens discovery (1=B06, 2=B08, 3= lTreu et all (|201lh . 4=this paper); Col. 13 is the 
band used for the lens modelling. 
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18.482175 


0.1153 


0.8710 


A 





,46 


4 


F435W/F814W 


J0825+2109 


126.320762 


21.161022 


0.1012 


0.3956 


X 





,26 


4 


F435W/F814W 


J0908+2730 


137.234813 


27.513829 


0.0219 


0.2831 


X 





,15 


4 


F160W 


J0915+4211 


138.81787 


42.19800 


0.078 


0.790 


A 





,59 


3 


F435W/F814W 


J0930+2855 


142.572560 


28.916732 


0.3051 


0.4594 


A 





,55 


4 


F435W/F814W 


J1000+2835 


150.011213 


28.595699 


0.0885 


0.8757 


X 





,41 


4 


F435W/F814W 


J1021 + 2028 


155.435084 


20.477728 


0.1208 


0.3507 


A 





,39 


4 


F435W/F814W 


Jllll+2234 


167.866760 


22.580729 


0.2223 


0.9887 


A 





,58 


4 


F435W/F814W 


J1135+3720 


173.77867 


37.33997 


0.162 


0.402 


A 





,32 


3 


F435W/F814W 


J1249+0225 


192.467683 


2.417972 


0.0965 


0.3517 


C 





,43 


4 


F435W/F814W 


J1403+1530 


210.799017 


15.511061 


0.0976 


1.0858 


C 






4 


F435W/F814W/F160W 


J1422+4134 


215.609046 


41.576612 


0.1011 


3.01 


B 





,20 


4 


F435W/F814W 


J1629+4708 


247.291334 


47.145309 


0.1282 


0.3739 


B 





,41 


4 


F435W/F814W 


J1633+1341 


248.414833 


13.687109 


0.1171 


0.2436 


C 





,31 


4 


F435W/F814W 


J2333-1042 


353.272047 


-10.703289 


0.1445 


0.1328 


B 





,37 


4 


F435W/F814W 



Table 2. Summary of basic properties of the SWELLS HST Cycle 18 targets. Col. 1 lists the lens ID; Cols. 2 and 3 the coordinates 
(J2000 in degrees); Cols. 4 and 5 give deflector and source redshifts; Col. 6 the lens grade (A=secure, B=probable ; C=p ossible, X=not a 
lens); Col. 7 the axis ratio as measured by SDSS; Col. 8 the reference for the discovery of the lens f3= lTreu et al.l (|201lh , 4=this paper); 
Col. 9 the HST filters. 



to n = 1/2 — i.e. we fit the lens galaxy with up to 10 
elliptical Gaussian components (this is simi l ar to the Multi- 
Gaussian Expansion method of ICappellaril d2002h ). For the 
more complex galaxies with star forming rings, we allow neg- 
ative Gaussians to be included as well. The star forming 
ring can be modelled as a smooth Gaussian-like profile with 
negative Gaussians suppressing interior surface brightness, 
resulting in a ring profile. We use the mult i- Gaussian model 
to subtract the lens galaxy from the image, resulting in an 



image of the lensed features only. These lens-subtracted im- 
ages are shown in the second column of Figure [H 



3.2 Gravitational Lens Model 

We fit Singular Isothermal Ellipsoid (SIE) lens models with 
external shear to the lensed features of each system. The 
source plane position (x s ,y 8 ) and the image plane position 
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ID 


B bulge 


-^bulge 


-^bulge 


-^bulge 


-^disc 


-^disc 


^disc 


-^disc 


-^e, bulge 
[arcsec] 


9bulge 


-^e,disc 
[arcsec] 


9disc 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 



SDSSJ0007+0053 


21 


79 


19 


12 








21 


79 


19 


72 











45 





68 


1 


40 





39 


SDSSJ0149-0010 


22 


96 


20 


09 








22 


93 


20 


56 











21 





56 





39 





48 


SDSSJ0329-0027 


22 


68 


19 


63 








21 


91 


19 


52 











16 





85 





90 





39 


SDSSJ0821 + 1025 


20 


69 


17 


57 








18 


75 


16 


81 








1 


28 





52 


2 


59 





15 


SDSSJ0822+1828 


19 


86 


17 


50 








20 


51 


18 


40 











65 





66 


1 


09 





24 


SDSSJ0825+2109 


48 


84 


20 


35 








18 


80 


16 


99 











32 





79 


1 


80 





24 


SDSSJ0915+4211 


19 


02 


16 


82 




15 


63 


19 


55 


17 


68 




16 


73 


1 


05 





89 


1 


40 





25 


SDSSJ0930+2855 


21 


53 


18 


28 








21 


97 


18 


95 











62 





53 


2 


36 





75 


SDSSJ1000+2835 


19 


56 


17 


32 








19 


61 


17 


71 











58 





79 


1 


43 





20 


SDSSJ1021+2028 


19 


75 


17 


42 








20 


82 


19 


13 











40 





44 





70 





17 


SDSSJ1111 + 2234 


20 


33 


17 


62 








19 


94 


17 


84 











71 





67 


1 


38 





46 


SDSSJ1135+3720 


21 


37 


18 


45 




17 


19 


19 


91 


17 


54 




16 


45 





56 





67 


1 


41 





26 


SDSSJ1249+0225 


20 


31 


18 


08 








19 


78 


17 


84 











37 





73 


1 


06 





32 


SDSSJ1403+1530 


20 


95 


18 


16 


16.88 






18 


98 


17 


02 


16.34 









60 





62 


1 


97 





36 


SDSSJ1422+4134 


48 


24 


19 


69 








18 


64 


16 


78 











54 





57 


2 


06 





16 


SDSSJ1629+4708 


20 


17 


17 


80 








19 


64 


17 


56 











39 





71 


1 


16 





25 


SDSSJ1633+1341 


21 


90 


19 


29 








19 


69 


17 


38 











91 





15 


2 


07 





26 


SDSSJ2333-1042 


20 


40 


16 


97 








19 


48 


17 


13 








1 


23 





65 


1 


32 





18 



Table 3. Structure of SWELLS HST Cycle 18 targets with multi-band imaging. Col. 1 lists the lens ID; Columns 2-9 the AB magnitudes 
of the bulge and the disc in the various bands, Col. 10 circularised half-light radius of the bulge; Col. 11 the axis ratio of the bulge; 
Col. 12 circularised half-light radius of the disc; Col. 13 the axis ratio of the disc. The stellar masses inferred from these magnitudes are 
listed in Table 0] 



ID logio( M *,bulge/M ) log lo (M^ disc /M ) log lo (M* ;bulge /M ) log 10 (M* ;disc /M ) 

Chabrier Chabrier Salpeter Salpeter 

(1) (2) (3) (4) (5) 



SDSSJ0007+0053 


10 


63 


± 





14 


10.26 


± 0.11 


10 


84 


± 





18 


10.43 


± 


0.13 


SDSSJ0149-0010 


10 


21 


± 





16 


9.90 


± 0.14 


10 


57 


± 





19 


10.21 


± 


0.13 


SDSSJ0329-0027 


10 


44 


± 





14 


10.26 


± 0.17 


10 


72 


± 





15 


10.57 


± 


0.15 


SDSSJ0821 + 1025 


10 


40 


± 





11 


10.60 


± 0.13 


10 


64 


± 





11 


10.78 


± 


0.16 


SDSSJ0822+1828 


10 


54 


± 





13 


10.09 


± 0.11 


10 


79 


± 





16 


10.34 


± 


0.16 


SDSSJ0825+2109 












10.58 


± 0.17 












10.81 


± 


0.16 


SDSSJ0915+4211 


10 


60 


± 





11 


10.17 


± 0.10 


10 


83 


± 





09 


10.43 


± 


0.09 


SDSSJ0930+2855 


11 


18 


± 





15 


10.85 


± 0.15 


11 


42 


± 





17 


11.11 


± 


0.17 


SDSSJ1000+2835 


10 


40 


± 





14 


10.18 


± 0.14 


10 


57 


± 





17 


10.34 


± 


0.14 


SDSSJ1021+2028 


10 


62 


± 





13 


9.74 


± 0.14 


10 


87 


± 





18 


9.98 


± 


0.12 


SDSSJ1111 + 2234 


11 


11 


± 





13 


10.88 


± 0.15 


11 


37 


± 





14 


11.10 


± 


0.15 


SDSSJ1135+3720 


10 


62 


± 





08 


10.91 


± 0.10 


10 


85 


± 





09 


11.09 


± 


0.07 


SDSSJ1249+0225 


10 


15 


± 





13 


10.16 


± 0.13 


10 


32 


± 





16 


10.50 


± 


0.15 


SDSSJ1403+1530 


10 


33 


± 





10 


10.51 


± 0.12 


10 


58 


± 





12 


10.76 


± 


0.10 


SDSSJ1422+4134 












10.56 


± 0.17 












10.84 


± 


0.15 


SDSSJ1629+4708 


10 


52 


± 





12 


10.49 


± 0.15 


10 


73 


± 





16 


10.78 


± 


0.12 


SDSSJ1633+1341 


9. 


88 




0. 


14 


10.67 


± 0.15 


10 


11 


± 





13 


10.84 


± 


0.20 


SDSSJ2333-1042 


11 


08 


± 





13 


10.86 


± 0.14 


11 


29 


± 





12 


11.07 


± 


0.15 



Table 4. Stellar mass estimates for the SWELLS HST Cycl e 18 ta rgets with multi-band imaging. Col. 1 lists the lens ID; Cols. 2 and 3 
t he stellar masse s of the bulge and disc assuming a Chabrier ( 2003) IMF; Cols. 4 and 5 the stellar masses of the bulge and disc assuming 
a lSalpeterl (119551 ) IMF. 



(x,y) of a ray are related by: 

x s = x - a x (x,y) 

y s =y-a.y{x,y) (1) 

where the deflection angles a x , y are given by 
iKeeton Kochanekl (|l998h . Denning <ip = y / q 2 x 2 + y 2 , 
where q is the axis ratio/flattening of the projected mass 
profile, the deflection angle formulae for the SIE component 
(i.e. before adding the external shear) are: 



ot x (x,y) 



a y (x,y) 



by/q 



: tan 



tanh 



(2) 



(3) 



as long as q < 1. If q > 1, then q can simply be replaced by 
q~ x and the orientation angle 9 rotated by 90 degrees, with 
the same result. We do not allow q = 1 but it can become 
arbitrarily close to 1 and therefore close-to-spherical lenses 
are not ruled out. The parameter b denotes the Einstein 
Radius of the lens, using the intermediate- axis convention. 

In many cases, we are attempting to fit this lens model 
to image data that do not contain very much information. 
Therefore, informative priors on the lens model parameters 
are required. The centre of the lens model is fixed to be 
within ± 3 pixels of the peak pixel of the light distribution. 
The position angle of the lens model is assigned an infor- 
mative prior, stating that the position angle of the lens is 
likely to be within ±10 degrees of the position angle of the 
lens galaxy light. The modelling results in an (intermedi- 
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ID 


(M/L) B 


(M/L) v 


(M/L) B 


(M/L) v 


(M/L) B 


(M/L) v 


(M/L) B 


(M/L) v 




Bulge 


Bulge 


Disc 


Disc 


Bulge 


Bulge 


Disc 


Disc 




Chabrier 


Chabrier 


Chabrier 


Chabrier 


Salpeter 


Salpeter 


Salpeter 


Salpeter 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


SDSSJ0007+0053 


2.37±0.70 


2.02±0.50 


1.37±0.33 


1.35±0.28 


3.98±1.54 


3.38±1.03 


2.05±0.55 


2.06±0.45 


SDSSJ0149-0010 


3.36±1.07 


2.65±0.71 


2.01±0.58 


1.76±0.40 


7.40±3.00 


5.58±1.82 


4.02±1.12 


3.55±0.78 


SDSSJ0329-0027 


3.38±0.95 


2.66±0.61 


1.65±0.65 


1.55±0.51 


6.62±2.06 


5.14±1.28 


3.28±1.02 


2.94±0.74 


SDSSJ0821 + 1025 


6.48±1.73 


4.52±1.05 


2.92±0.91 


2.44±0.64 


11.66±2.79 


8.20±1.67 


4.49±1.78 


3.82±1.18 


SDSSJ0822+1828 


3.45±1.14 


2.75±0.73 


2.48±0.68 


2.19±0.47 


6.30±2.28 


5.04±1.45 


4.49±1.70 


3.80±1.14 


SDSSJ0825+2109 






2.43±1.01 


2.11±0.67 






4.16±1.54 


3.62±1.05 


SDSSJ0915+4211 


5.16±1.40 


3.82±0.87 


3.54±0.97 


2.84±0.63 


8.45±2.11 


6.39±1.33 


6.41±1.58 


5.10±1.05 


SDSSJ0930+2855 


3.70±1.03 


2.84±0.65 


3.01±0.88 


2.37±0.55 


6.55±2.18 


4.98±1.35 


5.54±1.83 


4.38±1.17 


SDSSJ1000+2835 


4.12±1.33 


3.16±0.82 


2.98±0.99 


2.47±0.67 


6.33±2.69 


4.97±1.69 


4.36±1.49 


3.74±1.00 


SDSSJ1021+2028 


3.65±1.04 


2.85±0.64 


1.69±0.61 


1.60±0.45 


6.59±2.66 


5.19±1.72 


2.91±0.73 


2.77±0.59 


SDSSJ1111 + 2234 


3.23±0.85 


2.56±0.54 


1.87±0.60 


1.71±0.45 


5.97±2.00 


4.71±1.27 


3.08±0.88 


2.77±0.62 


SDSSJ1135+3720 


5.63±1.16 


4.05±0.73 


3.72±0.91 


2.92±0.62 


9.51±2.24 


6.86±1.41 


5.46±1.10 


4.41±0.73 


SDSSJ1249+0225 


3.74±1.11 


2.99±0.70 


2.53±0.85 


2.20±0.57 


5.49±2.08 


4.47±1.34 


5.65±1.92 


4.60±1.25 


SDSSJ1403+1530 


7.25±1.85 


5.00±1.12 


2.66±0.72 


2.25±0.48 


13.20±4.08 


9.06±2.43 


4.65±1.18 


3.90±0.83 


SDSSJ1422+4134 






2.22±0.92 


1.96±0.63 






4.28±1.55 


3.65±1.07 


SDSSJ1629+4708 


3.58±1.01 


2.83±0.64 


2.42±0.81 


2.06±0.56 


5.97±2.03 


4.72±1.33 


4.57±1.22 


3.88±0.85 


SDSSJ1633+1341 


4.17±1.24 


3.18±0.76 


3.83±1.21 


2.97±0.76 


7.14±2.10 


5.48±1.30 


5.80±2.41 


4.61±1.55 


SDSSJ2333-1042 


7.13±1.99 


4.92±1.21 


3.11±0.95 


2.51±0.61 


11.80±3.38 


8.20±2.03 


5.13±1.73 


4.21±1.14 



Table 5. Stellar mass-to-light ratios, in solar units (i.e. (M/M )/(L/L )), for the SWELLS HST Cycle 18 targets with multi-band 
imaging. Col. 1 lists the lens ID; Cols. 2-5 the mass-to-light ratio of the bulge and disc in the B and V bands assuming a IChabrierl (2003) 
IMF; Cols. 6-9 the mass-to-light ratio of the bulge and disc in the B and V bands assuming a lSalpeterl (|l955h IMF. 



SDSSJ091 5+4211, Kp 



Rotation Subtraction 




Figure 3. Rotation subtraction for the Kp image of the SWELLS 
system SDSSJ0915+4211. The lensed arc and counterimage are 
revealed by this process. The rotation-subtracted image is then 
used to mask the image so that a proper multi-Gaussian fit can 
be made. 



ate axis) Einstein radius b and an ellipticity parameter q for 
each system. Once the lens model parameters have been es- 
timated, the Einstein Radius in radians can be converted to 
an equivalent velocity dispersion asm using 



0"SIE 




(4) 



Ds 

47T " V ^ds 

wher e D s and Dds are the angular diameter distances (|Hogd 
Il999h from the observer to the source, and the deflector 
to the source, respectively. For simplicity, we modelled the 
sources with a single circular Gaussian component as this 
is adequate to explain the observed lensing configuration in 
most cases. For those cases where this model was inade- 
quate, the source was generalised to be either an elliptical 
Gaussian (e.g. for 0820) or multiple spherical Gaussians (e.g. 
for 0822). Due to the presence of systematics in the data 
from incomplete subtraction of lens galaxy light, and also 
the simplicity of the model used, we do not literally perform 
the inference with the lBrewer et al] (|2011al ) likelihood £(#), 
as this would result in inappropriately small uncertainties 
on the inferred quantities due to the systematic errors. Sim- 



ple, commonly adopted models for the likelihood function 
tend to assume that the measurement error in each pixel is 
independent, and that measurement error is the sole cause 
of discrepancies between the model and the data. This is not 
true in the presence of significant systematic effects such as 
residuals from an incomplete lens galaxy subtraction, and 
the assumption of a specific form (Gaussian) for the source 
light profile. 

To account for this and obtain realistic uncertainties on 
b and q despite the use of a simplified model, we adopted the 
following str ategy. This st r ategy has p reviously been briefly 
discussed bv lBrewer et al.l (|201ld ) and lBrewer et al.l (|2Q1 lal ) 
and will be presented in depth in a forthcoming contribution 
(Brewer et al 2012, in preparation). 

We carry out the sampling of the posterior distribu- 
tion for the lens model param e ters using a Nested Sam- 
pling algorithm (jSkillin el l200d : | Brewer, Partav, Csanvil 
l2011bh . Rather than simply sampling the posterior distribu- 
tion (which would yield unrealistic small uncertainties due 
to the presence of systematics and the simplified model) , we 
explore the dependence of the posterior distribution for the 
parameters on the temperature T: 



p(6\D-T) octt(0)C(0) 



l/T 



(5) 



Where denotes the parameters of the lens model, tt(0) is 
the prior and C(6) is the standard likelihood function. The 
temperature controls how conservative the posterior distri- 
bution is: T = oo reproduces the prior, and T — 1 gives the 
standard posterior obtained by ignoring the systematics. 

For each system, we then choose the highest tempera- 
ture T for which the posterior distribution over predicted 
lensed images (the 3rd panel in Figure [4|) matches the gross 
lensed features seen in the data. In other words, we force the 
posterior distribution to be as conservative as possible while 
still fitting the morphological features of the lensed images 
that are believed to be robust. Note that using a tempera- 
ture T is analogous to increasing the noise standard devia- 
tions on the data by a factor y/T. This procedure does not 
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change the values of our best estimates for the lens model 
parameters, but does increase the uncertainty around these 
estimates. Not carrying out this procedure would result in 
substantially underestimated uncertainties. 

See Figure [H for images of the lenses, along with galaxy 
subtracted images, lens models (chosen at random from the 
posterior distribution) and the residuals of the lens model 
fits. In several cases there are significant residuals remaining 
after the fitting procedure; these are generally due to the 
source having more structure than the surface brightness 
model that was assumed. In all cases the lensing aperture 
mass estimates are reliable. 



3.3 Notes on Individual Systems 

In this section we give brief notes on the image configura- 
tions of the lens models. 

0820+4847 Very little curvature is present in the arc. 
No definite counter-image is detected. An elliptical source is 
required, and the deflector is inferred to be very flat. 

0822+1828 Multiple (three) source components are re- 
quired. This results in tight constraints on the lens model 
parameters. 

0841+3824 Two source components are required. 
Whereas a bimodal posterior PDF is found, one mode chosen 
because of unrealistic small axis ratios for the other modes. 
A consistent lens model for this system has also been pre- 
sented by B08. 

0915+4211 Straightforward fit of the arc and counter- 
image with a single circular source. 

0930+2855 Significant residuals due to the assump- 
tions of no shear and isothermality (Auger et al., in prepa- 
ration) . The source has very high surface brightness and the 
system is a n example of an e arly- type /early- type lensing 
system fe.g. jAuger et al.ll201ll ). 

0955+0101 Straightforward fit of arc and counterim- 
age with a double image configuration. A consistent lens 
model for this system has also been presented by B08. 

1021+2028 We fit to the arc only. A bimodal solution 
was found, but we could rule out the mode which did predict 
a distant unseen counter- image. 

1029+0420 Classic double image configuration. A con- 
sistent lens model for this system has also been presented 
by B08. 

1032+5322 Arc and counterimage are well modelled 
by a single circular source. A consistent lens model for this 
system has also been presented by B08. 

1103+5322 Arc can be explained by a single circular 
source imaged four times. Lens model parameters are tightly 
constrained. A consistent lens model for this system has also 
been presented by B08. 

1111+2234 Arc can be explained by a single circular 
source imaged four times. Lens model parameters are tightly 
constrained. 

1117+4704 The single circular source is lensed into a 
quad configuration. Lens model parameters are tightly con- 
strained. 

1135+3720 The arc structure is complex, requiring 
two components for a good fit. Position of the counterimage 
is uncertain. Lens parameter constraints are not particularly 
tight. 



1203+2535 Position and morphology of the arc can 
be reproduced by a single circular source. Subtraction of 
the deflector galaxy was not straightforward for this system, 
leaving residuals that may add systematic uncertainties. 

1251-0208 Classic double image configuration. A con- 
sistent lens model for this system has also been presented 
by B08. 

1313+0506 Arc is well matched by a lensed circular 
source, but due to the non-detection of a counter- image, the 
constraints on the lens model parameters are weak. 

1331+3638 Classic cusp configuration leads to strong 
constraints on the lens model parameters. 

1422+4134 This is the only grade 'B' lens system for 
which a plausible model has currently been obtained. The 
two compact images are possibly double images of a back- 
ground QSO. This configuration yields weak constraints on 
the lens model, and the inferred stellar mass fractions place 
this system as a significant outlier from the relation in Fig- 
ure possibly due to its significant dust content or a po- 
tential misidentification of the source redshift. Hence, this 
system has been excluded from the IMF analysis of Sec- 
tion n 

1703+2451 Arc and counterimage can be explained by 
a single circular source. 

2141-0001 As discussed in paper II, the arc appears to 
be split into three images, but this is likely due to the pres- 
ence of two dust lanes. The arc morphology can be explained 
by a single circular source. The posterior distribution in- 
cludes models where the arc is split into distinct images, as 
well as models where the arc is continuous. 



4 COMPARING STELLAR MASS AND 
LENSING MASS 

In this section we derive lensing and stellar masses within 
the critical curves (§ I4.1[L and compare them to derive limits 
on the stellar in initial mass function (§ 14. 2p . 



4.1 Lensing and stellar masses within the critical 
curves 

As in Paper I, we fit the surface brightness models to the 
available high resolution imaging, i.e. the HST and Keck- AO 
imaging, using the multi-colour fitting code SPASMOID 
(Be nnert et al. [201 lh . The stellar mass (and surface bright- 
ness) model consists of an n — 4 Sersic profile (i.e. a de 
Vaucouleurs profile) for the bulge and an n — 1 exponential 
disc component. One advanta ge of SPASMOID over simi- 
lar methods such as GALFIT (IPeng et al.ll2002h include the 
ability to take into account prior information such as the fact 
that the disc component should have a larger effective ra- 
dius and a higher ellipticity than the bulge. Once we have 
obtained the surface photometry parameters, we then infer 
the stellar masses of the bulge and disc components using 
a Bayesian method for fitting stellar p opulation synthes is 
models to multi-band photometric data ([Auger et al.ll2009h . 
This method uses a grid of models output from a standard 
SPS code, and allows for marginalisation over unknown pa- 
rameters such as star formation history, metallicity, etc. A 
summary of the structural properties (stellar masses, sizes, 
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Figure 4. Images and gravitational lens models of SWELLS systems. Left to right: Image, lens subtracted image, lens model predicted 
image (for a model chosen at random from the posterior distribution), residuals of lens model fit. The black bar in the left panel has a 
length of 1 arc second. 
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Figure 4. (continued) Images and gravitational lens models of SWELLS systems. Left to right: Image, lens subtracted image, lens model 
predicted image (for a model chosen at random from the posterior distribution), residuals of lens model fit. The black bar in the left 
panel has a length of 1 arc second. 
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Figure 4. (continued) Images and gravitational lens models of SWELLS systems. Left to right: Image, lens subtracted image, lens model 
predicted image (for a model chosen at random from the posterior distribution), residuals of lens model fit. The black bar in the left 
panel has a length of 1 arc second. 
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Figure 4. (continued) Images and gravitational lens models of SWELLS systems. Left to right: Image, lens subtracted image, lens model 
predicted image (for a model chosen at random from the posterior distribution), residuals of lens model fit. The black bar in the left 
panel has a length of 1 arc second. 



axis ratios) of the SWELLS Cycle 18 targets with multi- 
band imaging is given in Tables [3] and |U For the other lenses 
we use the struct ural parameters as reported in SWELLS I 
(jTreu et al.ll201lh . 

We then integrate these models to obtain the total 
stellar mass M*,sps within the critical curve. This can be 
compared with the inferred mass within the critical curve 
from lensing (Mi ens ), yielding stellar mass fractions /* = 
M* 5 sps/Miens under both the Chabrier IMF assumption and 
the Salpeter IMF assumption. The reason for the choice of 
the critical curve as the integration aperture is that, heuris- 
tically, it is the aperture within which the m ass of the lens 
is most strongly constrained (e.g. , lTreull2010l , and references 
therein). To obtain uncertainties on the stellar mass frac- 
tions, we Monte- Carlo sample from the posterior distribu- 
tion for the lens model parameters and the posterior distri- 
butions for the total stellar mass. 

The resulting stellar mass fractions for the grade 'A' 
lenses, under the Chabrier and Salpeter assumptions, are 
plotted in Figure [5] as a function of lensing velocity disper- 
sion (a proxy for overall galaxy mass that can be measured 
accurately and independent of the size of the Einstein Ra- 
dius and the normalisation of the stellar IMF) . Remarkably, 
there is a smooth declining trend as a function of stellar 
velocity dispersion over more than a factor of two in stel- 
lar velocity dispersion (and a factor of ~10 in stellar mass). 
This trend can be due to two effects: i) a systematic vari- 
ation of the IMF; ii) a varying dark matter fraction within 
the Einstein Radius due to a change in stellar mass to dark 
matter ratio and the change in physical size of the Einstein 
radius. Disentangling the two effects is a major goal of the 
SWELLS survey and requires very realistic lensing and dy- 
namical models (e.g. paper II). However, as we will show in 
the rest of this paper, a very robust inference on the nor- 
malisation of the IMF at the low velocity dispersion end can 
be obtained based on a very simple argument. 

At the low velocity dispersion end, the Salpeter IMF 



assumption results in inferred stellar mass fractions greater 
than unity, which is clearly unphysical, especially consider- 
ing that there may be an additional contribution of cold gas 
distributed as the stars (paper II and references therein). 
This provides evidence against the Salpeter IMF and any 
IMF that would predict a similar or greater stellar mass. In 
Section [4. 21 we formalise this argument with a probabilistic 
model and derive quantitative constraints on the IMF for 
the lower mass (ctsie < 230 kms" 1 ) SWELLS galaxies. 



4.2 Constraints on the Initial Mass Function 

Our data do not allow us to probe the form of the IMF but 
they do allow us to place robust constraints on the integral 
of the IMF: if the assumed IMF yields a mass significantly 
greater than the inferred total mass from lensing (i.e., a stel- 
lar mass fraction much greater than unity) then this IMF 
is excluded by our data. We could assume that all of our 
galaxies have the same underlying central stellar mass frac- 
tion / t * rue , in which case the estimate for a Salpeter IMF (i.e., 
the variance- weighted mean of the points in Figure [5} would 
imply / t *ue > 1- However, there is certainly some amount 
of scatter in the central dark matter fraction and therefore 
in the stellar mass fraction. We therefore explicitly build a 
model for the true underlying stellar mass fraction that in- 
cludes Gaussian scatter and the fact that / t * rue is between 
and 1. In practice this is a truncated normal distribution, 
which we parameterise as having a peak at \i and charac- 
teristic width a. Then the probability of observing a set of 
galaxies with true stellar mass fractions {F*} is given by 



p({F*}\n,<T) = 11 -j 



exp[-^(F*- M ) 2 ] 



U ex p \-M F :-») 2 ] dF t 



(6) 



where F* = f* TU e,i £ [0> !]• We also introduce a parameter a 
that parameterises the difference between a true (unknown) 
normalisation of the IMF and the normalisation resulting 
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Figure 5. Stellar mass fraction within the critical curve vs lensing velocity dispersion for the SWELLS sample. The SLACS sample 
(with duplicates in SWELLS removed) ([Auger et al.ll2009h is also plotted for comparison. Note that, at the low mass end (to the left of 
the dashed line at <r SIE < 230 km s -1 ), the Salpeter IMF assumption results in inferred stellar mass fractions greater than unity, which 
is unphysical. 



from a Salpeter IMF: 

M*,true = Q:M*,Salp. 

This normalisation is always greater than 0, and a = 1 im- 
plies a Salpeter IMF (a Chabrier IMF typically would im- 
ply a ~ 0.55). At fixed total mass, an increase in a yields 
an increase in / t * rue , demonstrating that the parameters are 
highly degenerate. However, as noted previously, a cannot 
be raised to arbitrarily large values because the stellar mass 
fraction is bounded to be less than unity. In practice, a plays 
a similar role to a stellar mass-to-light ratio but has the ben- 
efit of being less dependent on dust and age due to having 
explicitly modelled the spectral energy distributions of the 
galaxies (see Section 

Given our data - total masses from lensing and stellar 
masses assuming a Salpeter-like normalisation of the IMF 
- and invoking Bayes' theorem allows us to write down the 
probability distribution for our unknown model parameters 
/x, cr, a and {F*}, given the data: 

p(a, fi, a, {F*}\{f*}) oc p(a,fj,,(T, {F*}) 

xp({/*}|a,/i,a,{^*}) (7) 
= p(ot)p(fi)p(a) 

x P ({F:}\^a)p({f:}\{F*}) (8) 

where /* = a:/* Salp are the scaled observed stellar mass 
fraction, {F*} is the set of true (unknown) stellar mass 
fractions, and the leading terms are priors on a, /i, and 
a. The term p({F*}\fi, a) is given by Equation \6\ and the 
term p({f*}\{F* }) is the likelihood of having observed the 



data, taken to be a multivariate normal distribution with 
means {F*} and standard deviations given by the measure- 
ment uncertainties (the galaxies are independent and there 
are therefore no covariances) . 

The inferences on the parameters (/x, a, a) were com- 
puted using Markov Chain Monte Carlo, with the true F*s 
explicitly integrated out: 



p(a,fi,a,\{f:}) = J p(a,^a,{F:}\{f:})d n F: (9) 

These inferences were calculated assuming a uniform 
prior on \i (a location parameter) between and 1, a log- 
uniform prior on a (a scale parameter) between 0.01 and 1, 
and a log-uniform prior on a (a scale parameter) between 0.2 
and 5, all independent. The results are shown in Figures [6] 
and Values of a less than about 0.9 are approximately 
equally plausible, but the evidence strongly disfavours values 
of a greater than unity. The evidence for the point hypoth- 
esis a = 0.5 is only about twice that of the point hypothesis 
a = 1, however the probability that the IMF is lighter than 
Salpeter is 98%. This result depends slightly on the chosen 
prior for /x, a and a, and becomes weaker if a greater prior 
probability is implicitly assigned to the hypothesis that the 
intrinsic F* values are clustered close to 1. For example, if 
the prior bounds on ji are extended to —0.1 and 1.1, allowing 
the true F* distribution to place most of its mass near 1, the 
posterior probability that the IMF is lighter than Salpeter 
is slightly reduced to 93%. The intrinsic scatter on the true 
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Figure 6. Inference on the IMF mismatch parameter a and the 
parameter the mean of the true distribution of stellar mass frac- 
tions. Values of a less than about 0.9 are approximately equally 
plausible, but the evidence strongly disfavours higher values of a. 
The vertical line denotes the Salpeter IMF, and the dashed line 
indicates a Chabrier-like (in terms of mass-to-light ratio) IMF. 
The probability that the IMF is lighter than Salpeter is 98%. 
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Figure 7. Inference on the parameter a, the intrinsic scatter 
in the true distribution of stellar mass fractions F* . The value 
of a is constrained to be low (< 0.16 with 95% probability), so 
there is evidence that the intrinsic stellar mass fractions are fairly 
homogeneous across the low mass end of the SWELLS sample. 



F* values is inferred to be small (a < 0.16 with 95% prob- 
ability) . 

Note that, if we take into account the fact that there 
can be an additional 20% of mass in cold gas distributed in 
the same way as the stars (as in Papers II), the probability 
that the IMF is lighter than Salpeter becomes greater than 
99%. Correspondingly, the 95% probability upper limit to a 
becomes 0.80. 



5 DISCUSSION AND COMPARISON WITH 
PREVIOUS WORK 

The main result of this paper is that the stellar mass of spiral 
galaxies with lens velocity dispersion between 160 and 230 
kms -1 (i.e. with masses comparable to or slightly above the 
Milky Way) is lighter than that predicted based on their 
colours using a Salpeter IMF. 

This conclusion is based on a robust and straightfor- 
ward argument, that is that the stellar mass cannot exceed 
the total mass measured by strong gravitational lensing. 
A high abundance of cold gas would make our result even 
stronger, requiring even less mass in stars and therefore a 
lighter normalisation of the IMF. Systematic uncertainties 
in the stellar mass estimates due to differences in stellar 
population synthesis models and a ssumptions are of the or- 
der of 0.05-0.1 dex for a fixed IMF ([Conrov. Gunn. fe White! 
120091 : lAuger et al.ll20Q9l : iTreu et alfeoioh and therefore not 
sufficient to modify our conclusions significantly. 

Even though our sample is the largest sample of spiral 
lens galaxies currently available it is important to consider 
whether this result could somehow be biased by unknown 
selection effects. In paper I we investigated the selection 
function in detail by comparing the size mass relation of 
SWELLS lenses to that of non-lens galaxies selected from 
SDSS and found them to be statistically equivalent. The 
detailed analysis of the lens SDSSJ2141-0001 presented in 
paper II indicate that that galaxy might be drawn from the 
densest part of the distribution of mass density profiles of 
spiral galaxies. Whether this is a generic feature of lensing 
selection needs to be verified by the detailed analysis of the 
entire sample. However, if lens galaxies were representative 
of the densest spiral galaxies, this would only make our in- 
ference stronger with respect to the overall population, since 
the densest galaxies are likely those with the highest stellar 
content. 

Our results are consistent with independent estimates 
of the normalisation of the stellar IMF of spiral galaxies 
based on a variety of techn iques. For example, ba sed on spi- 
ral galaxy rotation curves, iBell &; de Jona (|200ll ) find that 
a Salpeter IMF is too heavy and argue for a normalisation 
that is significantly lower. Based on the vertical velocity 
dispersion of discs, the Disk Mass survey (|Bershadv et al.l 
1 20 111 . 1 201 Oh argues that the mass-to-light ratio (and hence 
hence the normalisation of the stellar IMF) is significantly 
lower than that inferred assuming a Salpeter IMF. Inter- 
estingly, a qualitatively similar result is found for elliptical 
galaxies i n the same velocity dispersion range. Dynamical ar- 
gume nts (|Cappellari et al . 200 6; Dutton, Mendel, & Simard 
2012 ) as well as jo i nt dy namical and lensing arguments 
(jBarnabe et al.ll20lrtl201lh show that a Salpeter IMF is too 
heavy for low mass ellipticals. 

The light normalisation of the stellar IMF (or equiv- 
alently mass-to-light ratio) required for spirals and el- 
lipticals comparable to the mass of the Milky Way is 
in contrast with the heavier normalis ation required for 
more m assive systems. For exa m ple, [v an Dokkum & C onrovl 
d201Ch. ISpiniello et al.l (|201lh . lAuger et al.1 (|2Q10bh and 
ISonnenfeld et al.l (|201ll ). using different techniques and 
datasets, find that light Chabrier-like IMFs are ruled out for 
the most massive early- type galaxies, and a heavier IMF is 
preferred. Further evidence for this hypothesis is provided by 
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Treu et alJl20ld : lThomas et alJlioTH ) - and lNewman et al.1 
201lh for the case of massive brightest cluster galaxies - 
although those studies could not break completely the de- 
generacy between dark matter content and normalisation of 
the IMF. 

Taken together, these results are consistent with the 
hypothesis that the stellar IMF is not universal, but that it 
changes systematically with galaxy stellar velocity disper- 
sion and mass (and therefore halo mass to first approxima- 
tion). Since stellar velocity dispersion correlates with the age 
of the stellar populations, this non-universality might reflect 
the non-universality of the conditions that were prevalent 
when the stars were formed. The obvious suspects are the 
density of the gas, and the condition of ambient background 
radiation. A key open question is whether this variation in 
the IMF can be understood theoretically based on first prin- 
ciples. 

It is important to emphasise that lensing and dynam- 
ical measurements constrain the total normalisation of the 
IMF and not its specific form. Therefore, by themselves they 
cannot constrain the shape of the IMF and reveal whether 
the change in the normalisation is due to a change in the 
abundance of low mass stars or remnants of high mass stars 
or any combination thereof. Independent constraints such as 
those from stellar populations diagnostics or from galaxies 
at different redshifts, probing different ranges in stellar mass 
are necessary to break this degeneracy in the interpretation 
of the trend. 

A final issue that we have not addressed in this paper is 
how much of the trend shown in Figure [5] is due to a varying 
IMF and how much is it due to varying dark matter con- 
tent. For the higher velocity dispersion SLACS galaxies weak 
and str ong lensing and dynamical data allowed I Auger et al.l 
(2010b) to disentangle the two effects, showing that the most 
likely explanation is that both the dark matter content and 
the norm alisation of the IMF increase with galaxy mass 
fsee also | Cardone et aD 120091: ITortora et all |2Q10L l201ll: 
Deason et al.1 1201 ll: lGrillo"fe Gobat 201 (TTrh omas et al.l 



201 ll : button et aTfeoilal : ITortora et aI1l2Q12h . 

Future papers of this series (Barnabe et al. 2012; Dut- 
ton et al. 2012, in preparation) will combine the lensing 
information with dynamical models of kinematic observa- 
tions to construct two component mass models and disen- 
tangle the luminous and dark components for the SWELLS 
sample as well (see button et al.1 2011 bl : iTrott et al.ll2010l : 
Ivan de Yen et al.ll2010l : ISuvu et al.ll201lh . This will allow us 
to differentiate trends in the stellar IMF from trends in dark 
matter content over the entire range in velocity dispersion. 



6 SUMMARY AND CONCLUSIONS 

In this paper we have presented complete multiband HST 
imaging (including new images obtained as part of our cycle 
18 HST program, GO-12292) for the largest currently known 
sample of gravitational lens systems where the deflector is 
an inclined discy galaxy. The SWELLS sample consists to 
date of 20 secure lenses and 6 probable lenses. 

For each secure lens and one probable lens, we have 
derived gravitational lens models based on a single isother- 
mal ellipsoid deflector plus external shear. In addition we 
have used the multicolour HST imaging to estimate the stel- 



lar mass of the bulge and disc component of each deflector 
based on stellar population synthesis models. We have used 
the two measurements of stellar and total mass within the 
critical curve to study the normalisation of the stellar IMF 
for spiral lens galaxies. 

Our main conclusions can be summarised as follows. 

• The ratio between stellar mass and total mass within 
the critical curve correlates with the velocity dispersion of 
the deflector galaxy. By combining the SWELLS sample 
with the SLACS sample of more massive early-type galaxies 
we find that the trend extends between 160 and 350 kins" 1 . 
This trend can be due to variations in the normalisation of 
the stellar IMF and/or in non baryonic dark matter content. 

• At the low velocity dispersion end (below 230 kms -1 ) 
the Salpeter IMF predicts, unphysically, that the mass in 
stars is greater than the total mass. 

• Based on a rigorous probabilistic model, we find that 
the probability that the true stellar masses are lighter than 
implied by the Salpeter IMF is 98% for galaxies with velocity 
dispersion below 230 kms -1 , neglecting the contribution of 
cold gas. If the contribution of cold gas is taken into account 
the probability increases to greater than 99%, and a 95% 
upper limit on the stellar masses is 80% of the mass implied 
by the Salpeter IMF. 

• The light IMF required for low mass spiral galaxies is at 
variance with the Salpeter or heavier IMF required for more 
massive early- type galaxies, consistent with a non universal 
stellar IMF. 
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APPENDIX A: COMPARISON OF LENS 
MODELS WITH PREVIOUS MODELS 

The SWELLS sample contains six spiral lens galaxies that 
were originally discovered and modelled in SLACS. As a 
sanity check, we compare the inferred Einstein Radii of these 
lenses from the models p r esent ed in th is paper with those 
reported by [AugerelaD (|2009h . Since lAuger et al.l (|2009h 
reported Einstein Radii in units of kiloparsecs, we convert 
our Einstein Radii to kiloparsecs also. The comparison is 
shown in Figure IA1I and demonstrates agreement between 
the two sets of models. 
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